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The damping of single-particle and collective motion in exotic isotopes is a new topic 
and its study may shed light on basic problems of nuclear dynamics. For instance, it is 
known that nuclear structure calculations are not able, as a rule, to account completely 
^ ■ for the empirical single-particle damping. In this contribution, we present calculations of 
the single-particle self-energy in the case of the neutron- rich light nucleus 28 0, by taking 
proper care of the continuum, and we show that there are important differences with the 
case of nuclei along the valley of stability. 

1. Introduction 
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A useful approximation to describe atomic nuclei is in terms of an average potential 
in which neutrons and protons move independently of each other. This average potential 
is not a static one, and undergoes conspicuous fluctuations which give rise to collective 
oscillations of the system. The coupling of the particle motion to the collective modes is 
at the origin of the energy dependence of the average potential (cf., e.g., Refs. [ |T], |2[). The 
empirical energy-dependent average potential is also able to account (at positive energy) 
for the behavior of a particle impinging on the nucleus in a scattering experiment, as 
it includes an imaginary part which corresponds to the absorptive scattering channels. 
Within the so-called nuclear structure approach, the Hartree-Fock (HF) mean field can 
thus be viewed as a starting point to which corrections can be added, the leading ones 
being associated with processes in which a nucleon excites 1 particle-1 hole (lp-lh) pairs 
- preferably correlated pairs, that is, nuclear collective vibrations. Such corrections have 
been extensively studied for nuclei along the valley of stability. As a rule one finds that 
this approach is able to account only for about 50% of the empirical damping of single- 
particle motion [ Q . This is still an unresolved issue and points to the limits of accuracy 
of our understanding of the nuclear dynamics. 

In keeping with these facts, it is quite natural to ask the question whether in the recently 
discovered exotic nuclei (for instance, in systems with larger neutron excess compared to 
the stable isotopes) the situation described above is changed or not. This is one motivation 
for the present work. An even stronger motivation comes from the fact that, as it is well 
known [ [3| , the coupling of the single-particle motion to nuclear vibrations is responsible 
for most of the spreading widths of giant resonances in stable nuclei. In the case of exotic 
nuclei, the importance of the coupling with the continuum, that is, of the escape width 
of both the single-particle states and collective states has been stressed in many works. 
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On the other hand, not much work has been done to calculate the spreading width due 
to coupling with more complicated configurations (2p-lh for single-particle states, 2p-2h 
for giant resonances). 

In a previous work [ |J , we have calculated the quadrupole response in the neutron- rich 
nucleus 28 by taking into account both the escape and spreading width. We have reached 
the conclusion that the spreading width is as important as the escape width, at least in 
the example considered. Some of the results of Ref. [ |J were considered to be model- 
dependent, so that a more refined calculation is highly desirable. In fact, so far one has 
been able to take properly into account the continuum part of the single-particle spectrum 
in calculations at the lp-lh level and "more refined calculation" means therefore, within 
the present context, a proper treatment of the continuum also at the 2p-2h level. With 
this purpose we performed calculations in which the self-energy of particles (holes) is 
obtained by coupling with 2p-lh (2h-lp) configurations and for the first time in the study 
of exotic nuclei, we properly treat the continuum at this 2p-lh (2h-lp) level. 

Preliminary results are reported in this contribution. In a previous work [ H], con- 
siderations based on analytical expressions had been presented and we confirm here the 
soundness of these considerations by means of some numerical findings. 

2. Imaginary part of the single-particle self-energy in exotic nuclei 

The imaginary part of the single-particle self-energy felt by a nucleon in the presence 
of an A-particle bound system has been one of the main concerns of many studies in 
nuclear structure, both experimental and theoretical. When the nucleon is at positive 
energy and its wave function has ingoing wave boundary conditions, this imaginary part 
corresponds to the imaginary part of the optical potential. Whereas for heavy ions the 
dominant contribution to the optical potential comes from transfer phenomena and this 
can be accurately calculated [ in the case of nucleon scattering the theoretical under- 
standing of the optical potential is still a rather open problem. In fact, there is not a 
clear-cut explanation why, for instance, the optical potential calculated by starting from 
nuclear matter studies and using a local density approximation to deal with finite nuclei, 
reproduces more or less the strength of the empirical optical potential; and why, on the 
other hand, within the above mentioned nuclear structure approach, only about one half 
of the potential strength can be accounted for. 

Within this latter approach [ 0, H], the optical potential is calculated by adding to the 
HF mean field (Vhf) the contribution AV coming from the excitation of 2p-lh states 
(mainly correlated, that is, lp-1 collective vibration). Only AV gives rise to an imaginary 
part in the single-particle potential. In particular, only the upper two diagrams associated 
with AV and shown in Fig. 1, are relevant. We report the analytic expression for the 
upper left diagram, contributing to the particle self-energy, which reads 

r tt ( , w r WA) (r) u^-7 Wx \r') 
ImAVij(r,r';uj) — / duj\ (—ir)v(r)v(r') — — x 

i',j',>> r r 

x ^)(r)^)(rO (Zil ^y )2 . (1) 

The quantum numbers l,j (l',j f ) refer to the initial (intermediate) state, u(r) are radial 
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Figure 1. The diagrams corresponding to the leading terms of the particle (left column) 
or hole (right column) self-energy. 



wave functions of the particles, Sg(r) are transition densities of lp-lh pairs coupled to 
multipolarity A and v (r) is the p-h interaction derived from the effective force from which 
the HF mean field is determined (in our case, we have used the Skyrme force SIII). The 
link between Eq. (|ip and the empirical imaginary part of the optical potential W is as 
follows, 

W(r, r'; uj) = V -^l/mAVj,- (r, r'; w) . (2) 
; . An 

The real part is obtained in a similar way, or through dispersion relation techniques 
(see Ref. [ [I], |2|] and references therein). A procedure exists to connect the non-local 
potential with its local equivalent, which is the quantity to be compared with the empirical 
parametrizations or to be given as an input in a reaction calculation, e.g, of DWBA type. 

The novel feature of Eq. (|l|), in comparison with what has been used in the past, is 
that we treat properly the continuum also at the level of the 2p-lh doorway states. In the 
calculations performed in the seventies in well-bound nuclei like 208 Pb, discrete particle 
states were employed as the whole system was set in a box and an averaging parameter was 
employed in Eq. ([!]) to ensure the match of the initial lp and intermediate 2p-lh energies 
(see Ref. [ |IJ and references therein). This procedure is satisfactory for well-bound systems 
as we have been able to confirm by making use of Eq. (^ for a test calculation in 208 Pb 
where we have reproduced essentially the results of Ref. [ ||. For systems with loosely 
bound nucleons, where the correct treatment of the continuum is more important, the 
use of discrete particle states is not able to reproduce the results that we are going to 
illustrate in Sec. [| 
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Before that, we use Eq. ([[]) to get some insight in the qualitative low-energy behavior 
of the imaginary part of AV (by "low-energy" we mean for values of the energy uj close to 
the particle emission threshold uJth)- We use the fact that in the case of nuclei with loosely 
bound neutrons, the low-energy part of the (lp-lh) multipole response is characterized 
by a pronounced "threshold effect" , that is, by a sudden increase of the strength function 
S{uj) above uj th [ [| H. On ground of simple arguments, it is expected that 

SH~( W -%) ! ' +1 / 2 , (3) 

where I' is the orbital angular momentum of the particle states contributing to S(uj). 

Let us consider a single term of the sum appearing in eq. flU) and fix r = r'. A plane- 
wave approximation for tt^r WA ' ) (r) suggests that it contributes with a factor k 1 ^^ 2 if 

kpart = hT l ^J2m{uj — oj\) is close to zero. Adding the condition on the normalization of 
the radial transition densities given by 

| J dr r 2+x 5g M (r)\ 2 = S(uj), 

with the asymptotic behaviour of the strength recalled in Eq. ([3]), we obtain 

ImAV ~ / duj x (uj - uj x ) v+1 ' 2 {uj x - uJ th ) l+l ' 2 . (4) 

If uj is close to uj t h, i.e., uj = uj t h + 8 and uj\ = uj t h + 8/2, we find that the imaginary 
part of the optical potential behaves approximately like 5 l +l +2 , that is, a fast increase 
just after threshold which has no counterpart in stable nuclei (where ImAV ~ {uj — ep) n , 
with 1 < n < 2 [ 0) and has consequences on the qualitative features of scattering 
experiments with exotic beams. A realistically calculated ImAV follows indeed this 
asymptotic behaviour and this is shown in the next Section. 

3. Results 

We have chosen as an example of neutron-rich nucleus, the isotope 28 O. It is double 
magic and the neutron separation energy (about 1 MeV according to the HF calculation 
already mentioned, employing the SIII interaction) lies between the stable nuclei and 
the extreme cases of nuclei with a very extended halo like 11 Li. We must recall that a 
number of experimental evidences have been recently accumulated which point to the non- 
existence of this isotope as a bound system [ [TOfl . Since essentially all calculations predict 
28 to be bound, a strong point should be made about the necessity of improving our 
tools for nuclear structure studies, to face the outcome of the considerable experimental 
efforts made in the field of drip-line nuclei. 

We have calculated, in the case of neutrons, the quantity W defined by Eq. (|2]). First, 
we have solved the HF equations with the Skyrme interaction SIII on a radial mesh of 
250 points with a 0.1 fm step. As a result, the levels shown in Fig. 1 of Ref. [ ^] are 
obtained. In particular, the highest occupied neutron orbital d.3/2 has a small binding 
energy of about 1.1 MeV while the protons are all bound by more than 30 MeV. This 
large asymmetry between the neutron and proton mean fields is clearly due to the large 
neutron excess. Its consequences on the multipole response are quite dramatic as remarked 
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Figure 2. Strength function for the isoscalar 2 + , 3 and isovector 1 , calculated within 
continuum- RPA. 



first in Ref. [ [J. In fact, we have calculated [ [|, |5] the strength function associated with 
isoscalar quadrupole, octupole and dipole operators, within the framework of continuum 
RPA in coordinate space [ PH . This theory is able to treat properly the particle states 
in the continuum, because the one-body propagator has an exactly computable form 
if expressed as function of r, r '. We have kept 25 fm as the upper limit of the radial 
integrals while the radial mesh has a 0.5 fm step. The strength functions mentioned above 
are depicted in Fig. [2|. Their main feature is that they reflect how "excess neutrons" and 
"core particles" degrees of freedom are decoupled in systems like the one at hand. The 
low-energy part, although characterized by a strong enhancement just above the particle 
threshold, is not associated to any collective effect. What is enhanced so much to form a 
well-defined bump, are the single-particle transitions from the d 3 / 2 neutron state to s-, p-, 
or d-states in the continuum: because of the small d3/2 binding energy, the wave function 
of these neutrons is so extended that its overlap with continuum wave functions is large 
and this makes the matrix elements of the multipole operators quite large. On the other 
hand, the higher energy region of the strength functions is characterized by states to 
which more ph configurations partecipate (neutron excess states play a predominant role 
but neutron core states are not negligible). Transitions from proton states are completely 
decoupled and lie above 30 MeV. 

From the continuum RPA calculation we can extract the radial transition densities 
Sg^ x \r) which appear in Eq. ([]]). For each multipolarity A = 2 + , 3~, l - , we take energy 
bins of 1 MeV for uo\ and we associate to each of them the appropriate Sg(r) so that 
the integral / duj\ can be performed. The continuum wave functions u(r) are obtained 
by integrating the radial Schrodinger equation including the kinetic energy plus the HF 
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Figure 3. Optical potential W as defined in Eq. (H), plotted in the case of r = r' — 3.2 
fm, as a function of the energy uj close to u t h (thick line). In this low-energy region, the 
approximation W ~ (uj — Uth) 4 is valid (thin line). 



potential at positive energy uj — u>\. 

We first study the low-energy asymptotic behavior of W (Eq. (H)), to see whether it 
meets the expectations from Eq. (^). In Fig. ^ we show (by means of thick line) the values 
of W for r = r' fixed approximately at the nuclear surfaceQ The thin lines is a function 
of the type C(u — uj t h) 4 , with C constant and chosen to get the best fit to the microscopic 
calculation. The good agreement between the two lines confirms the soundness of the 
argument of the last Section, about the asymptotic behavior of the optical potential. 

In Fig. |], W(r,r f = r;u) is plotted as a function of r for three different values of the 
energy uj. As in stable nuclei, this function is surface peaked. On the other hand, the 
width of the peak is rather large, as can be expected if the nucleus has extended radial 
wave functions. All these characteristics of the imaginary part of the self-energy in exotic 
nuclei will be studied in more details in a forthcoming work. We can also notice that 
for uj = 25 MeV a second bump, at the interior of the nucleus, shows up. This can be 
interpreted as follows. In the calculation of the terms of the sum appearing in Eq. (|2|), one 
performs integrals of the type f£ duj\. If uj is large, one includes, through the transition 
densities, contributions from protons and these are of course peaked elsewhere than on 
the nuclear surface, since the protons are concentrated in a much smaller region. In this 
respect, this interior peak is another example of the general statement according to which 
core particles and excess particles are decoupled in this light, exotic isotopes. 

1 Fig. 1 of Ref. [ |j| shows that the presence of nucleons with binding energies of a few MeV give rise to a 
very diffuse surface. We have taken here a value of 3.2 fm, where the density is approximately half of its 
central value. 
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Figure 4. Optical potential W as defined in Eq. (|2|), plotted in the case of r = r', as a 
function of the radial distance for three different values of the energy uj (dotted line: 5 
MeV; dashed line: 15 MeV; full line: 25 MeV). 



4. Real part of the self-energy 

The off-shell diagrams which are shown in the lower part of Fig. [I] contribute to the 
real part of the self-energy, in addition to the two which appear above. It has been shown 
in the past (see Ref. [ |I[ and references therein) that dispersion relation techniques can 
be used to determine the values of Re AV, once Im AV has been calculated. We do not 
recall here but the basic idea, by using again the particle self-energy (diagrams in the left 
column of Fig. |I]). 

The imaginary part of the self-energy (Eq. (p])) can be obtained as a limit, 
Im/Wi;{r,r';uj) = Jmlim / du\ [ duj' ; — r , (5) 

where / = fij-,i f ,j',x(r, r'; uj, uj', u\) denotes the product of the two particle-vibration ver- 
tices. The integral over duj\ runs from uj t h to +oo. The integral over duj' can be extended 
from — oo to +oo and then the contribution from the 2 particles-1 hole-1 phonon interme- 
diate states (lower left diagram of Fig. H) is also included. This is just a formal extension 
if one takes the imaginary part of the limit rj — > 0, since this results in a delta-function 
8{u' — {uj which rules out the 2 particles-1 hole-1 phonon states contribution if to is 

positive. But if we consider the real part in the limit rj — > 0, we easily get the expression 
for the real part of the self-energy, 

ReAVij(r,r';uj) = Re lim / duj\ ( duo 1 ' 



uj - (u 1 + lux) + irj 
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, ImAVij(r, r'; u') 
uj — UJ 1 



(6) 



In the case of 28 0, by using ImAVij(r, r'; u') calculated as it is described in Sec. 3, we 
have obtained results also for the real part of the self-energy given by Eq. @. The problem 
we have addressed is actually the energy shift due to the coupling with the states labelled 
by A, U\ (which are collective vibrations if uj\ is around 15 MeV and are single-particle 
transitions at lower energy). This shift is given exactly by 

AE nlj = J drdr' ru n ij{r) ReAVij(r,r';uj) r'u n ij(r'), (7) 

where ui is usually fixed as the unperturbed (i.e., HF) energy E^- of the level under 
consideration. The expressions for ReAVij in the case of hole states are analogous to 
those derived above for the case of particles. In the case of 28 O calculated within the 
SIII-HF procedure, no unoccupied particle states can be found at negative energy and 
the only particle resonance at positive energy is the £7/2, so that the particle spectrum 
consists essentially only of a smooth continuum. Therefore, we have calculated the shift 
AE n ij for neutron hole states, in particular for the loosely bound d.3/2 orbital. 

We find a positive shift of 330 keV, to be compared with E^ = -1.1 MeV. This makes 
the nucleus more close to being unbound and of course would in turn affect the strength 
function of the multipoles we have considered in Sec. 3. As a conclusion to this part, the 
following two statements about the role of the real part of the self-energy in neutron-rich 
nuclei can be made, (a) The particle shift of 330 keV obtained for the loosely bound 
d3/2 orbital forces us to ask questions about the reliability of simple mean field methods 
to predict the single-particle energies and consequently the position of drip lines, as well 
as the properties of the excited states. This is not the case for stable nuclei, since one 
usually obtains corrections of the order of less than 1 MeV to the energy of states which 
are usually 7-10 MeV bound. These corrections are important to give the correct density 
of states around the Fermi energy and effective mass but do not alter the predictive power 
of theories like the standard RPA in which these corrections are not taken into account, 
(b) The fact that the net result of 330 keV comes from a strong cancellation between the 
two diagrams of the right column of Fig. [I] (930 keV - 600 keV), one has to wonder about 
the use of second RPA (as we did in [ ^]) where only the first of the two diagrammatic 
contributions is included. In conclusion, we believe that further research along this line 
can address some basic problems of the present physics of nuclei with large neutron excess. 
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